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ABSTRACT 

Context. The Rossby wave instability in astrophysical disks is as a potentially important mechanism for driving angular 
momentum transport in disks. 

Aims. We aim to understand this instability in an approximate three-dimensional disk model environment which we 
assume to be a single homentropic annular layer we analyze using disk shallow-water theory. 

Methods. We consider the normal mode stability analysis of two kinds of radial profiles of the mean potential vorticity: 
The first type is a single step and the second kind is a symmetrical step of finite width describing either a localized 
depression or peak of the mean potential vorticity. 

Results. For single potential vorticity steps we find there is no instability. There is no instability when the symmetric 
step is a localized peak. However, the Rossby wave instability occurs when the symmetrical step profile is a depression, 
which, in turn, corresponds to localized peaks in the mean enthalpy profile. This is in qualitative agreement with 
previous two-dimensional investigations of the instability. For all potential vorticity depressions, instability occurs for 
regions narrower than some maximum radial length scale. We interpret the instability as resulting from the interaction 
of at least two Rossby edgewaves. 

Conclusions. We identify the Rossby wave instability in the restricted three-dimensional framework of disk shallow water 
theory. Additional examinations of generalized barotropic flows are needed. Viewing disk vortical instabilities from the 
conceptual perspective of interacting edgewaves can be useful. 

Key words. Hydrodynamics, Astrophysical Disks - theory, instabilities 



1. Introduction 

The Rossby wave instability (hereafter RWI) is a promising 
candidate mechanism to account for the observed anoma- 
lous transport of cold astrophysical disks Q. The process in- 
volves the instability of waves in an environment in which 
there are radial variations in a potential vorticity quantity 
(PV hereafter) . This effect is the disk analog of the Rayleigh 
instability of stably stratified barotropic shear flows that 
are familiar in meteorology and geophysical fluid dynam- 
ics. The original works wherein the RWI was proposed for 
disks (Lovelace et al. 1999; Li et al. 2000; Li et al. 2001) 
demonstrated the existence and evolution of the instabil- 
ity in vertically integrated disk models. The earliest known 
(to this author) identification of the relevance of PV dy- 
namics to astrophysical settings is found in Lovelace & 
Hohlfeld (1978) containing a study of the stability of cold 
self-gravitating galactic disks subjected to radially local- 
ized perturbations. For the galactic disk study as well as 
the RWI, instability requires the existence of a minimum 
or maximum in the associated radial PV-profile of the disk0 
Because of the potential importance of this mechanism 
for disks, it is worthwhile to consider and understand this 
instability in disk settings, which are three dimensional, at 
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1 i.e., non-magnetized disks. 

2 Note that in Lovelace & Hohlfeld (1978) the associated quan- 
tity f(r) corresponds roughly to the inverse of the PV. 



least in part. The goal of this study is to expose the ma- 
chinery of this process as clearly as possible in the context 
of a three-dimensional theory. 

Disk shallow water theory (Umurhan 2008, and here- 
after DSW-theory) is a model reduction of the disk equa- 
tions that is three-dimensional but asymptotic in the sense 
that azimuthal scales are much larger than the correspond- 
ing vertical and radial scales. It is essentially a model en- 
vironment representing vortex dynamics occurring on thin 
annular sections of disks over timescales which are much 
longer than the local disk rotation time. The original mo- 
tivation for this approximation was to develop a three- 
dimensional framework to describe the dynamics of elon- 
gated vortices known to emerge in two-dimensional studies 
such as that reported in Godon & Livio (1999). The bene- 
fit of using the DSW-theory is that it allows a transparent 
analysis of vortex dynamics free of other physical processes 
that are likely not to play a significant role - in particular, 
free of both acoustic and gravity mode oscillations. As such, 
by means of this framework one may develop a mechanical 
understanding of disk-related vortex dynamics. 

In Umurhan (2008), DSW-theory was applied to under- 
stand the emergence of one form of the strato-rotational 
instability, which is an instability of a stably stratified 
Rayleigh stable shear profile in a channel (Yavneh 2001; 
Dubrulle et al. 2005). The no-normal flow boundary condi- 
tions at the walls of the domain bring into existence edge- 
waves (Goldreich et al. 1986) that propagate along the walls 
in opposite directions with respect to one another. However, 
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because the edgewaves may also interact with one another 
across the domain, if the separation of the walls fall within 
a specific allowed range, then the phase speeds of the edge- 
waves can 'resonate' with each other and become unstable 
(Umurhan 2006; Umurhan 2008). 

This edgewave dynamical picture has been used to un- 
derstand the development of many forms of barotropic geo- 
physical fluid instabilities including, among others, the orig- 
inal problem considered by Rayleigh (Hoskins et al. 1985; 
Baines & Mitsudera 1994; Heifetz et al. 1999) as well as 
astrophysical processes such as the Papaloizou-Pringle in- 
stability (Papaloizou & Pringle 1984; Goldreich et al. 1985). 
In the geophysical fluid dynamics context, this is frequently 
referred to as the instability of counter-propagating Rossby 
waves (e.g., Heifetz et al. 1999). The use of edgewaves as an 
interpretive tool has precedence in the analysis of plasma 
instabilities. For example, the diocotron effect, which is 
an instability associated with finite thickness charge sheets 
subjected to E x B drift, is rationalized to arise from the 
interaction of edgewaves counter-propagating along the sur- 
faces of the charge sheet (Knauer 1966). 

DSW-theory is an appropriate platform to examine 
the RWI. Indeed, an examination of the vortex structures 
emerging from nonlinear calculations of the instability (see 
Fig. 1 of Li et al. 2000, hereafter LFLC-2000) shows that 
the coherent vortex structure is severely elongated in the 
azimuthal direction while being radially tightly confined. 
This structural quality is similar to the vortices reported 
in Godon and Livio (1999). The advantage of DSW-theory 
is that it can capture some essence of three-dimensionality, 
while its disadvantage is that, because of its construction, 
small azimuthal scales cannot be represented. By contrast, 
the vertically integrated model environment used in the 
original RWI studies are able to capture small azimuthal 
scales while coming at the expense of representing three- 
dimensionality. 

In this study, the RWI is examined in the framework 
of DSW-theory. The stability of two analytically tractable 
mean PV-profiles is studied: (i) a single defect, or 'single- 
step', and (ii) a double defect, or 'symmetric-step' (see Fig. 
[1]). As in the case of the stratorotational instability, the 
RWI is shown to result from the interaction of Rossby edge- 
waves that propagate along the locations of the steps of the 
mean PV-profile. However, the instability manifests itself 
only in the symmetric step case. In Section 2, the equa- 
tions of motion and their perturbations are presented. In 
Section 3, the equations are analyzed for two kinds of mean 
potential vorticity profiles including the one just described. 
Analytic solutions are determined throughout by consider- 
ing the quasi-linear problem posed by considering the spe- 
cial case 7 = 3, where 7 is the ratio of specific heats. Section 
4 concludes with comments and a discussion about recent 
three-dimensional RWI calculations (Meheut et al. 2010). 

2. Equations and potential vorticity conservation 

The asymptotic equations upon which DSW-theory is 
based, describing the evolution of annular disk sections, 
was developed in Umurhan (2008). In that study, the equa- 
tions of motion for an annular section in which the ver- 
tical/radial extent C is much smaller that the azimuthal 
scale R was developed as an asymptotic series expansion in 
powers of e = C/R. The timescales of the dynamics (T) is 
comparatively long with respect to the local rotation time 
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Fig. 1. Sketches of profile plots considered in this work, 
(a) Single defect/jump located at x = 0. (b) Double defect 
(symmetric step) with jumps located at x = ±A. Double 
defect depicted is a depression whose mean (midplane) en- 
thalpy profile corresponds to a peak. 

(1/ O ) of the annular disk centered on R i.e., T ~ l/(eflo). 
To ensure a non-trivial asymptotic balance, the radial and 
vertical velocities viewed in the reference frame of the an- 
nulus rotating with fig must scale as ec s where c s is the 
sound speed of the gas, while the corresponding deviations 
of the azimuthal velocities must scale as c s . This dispar- 
ity in aspect and velocity ratios is an asymptotic descrip- 
tion of azimuthally elongated slowly overturning vortices. 
Furthermore, the reduced equations of motion in the an- 
nular disk section possess quasi-geostrophic characteristics 
familiar in meteorological studies. 

Single constant entropy layers are considered when the 
equation of state is polytropic, i.e., it is assumed that 
P = Kp 1 everywhere at all times. This corresponds to a 
situation in which the specific entropy inside, defined by 
S = C v \n(P/p~<) — C v In if, is everywhere a constant. This 
homentropic layer is considered in the following. If there are 
radial variations of the mean height of the disk (h) , then the 
vertically integrated entropy of the disk also exhibits radial 
variations even though the specific entropy is constant. This 
point is made to keep in mind how the equations analyzed 
here compares to the vertically integrated two-dimensional 
equations studied in Li et al. (2000). 

After non-dimensionalizing all the quantities according 
to the above, the leading order equations describing the evo- 
lution of varicose disturbances of a single constant specific 
entropy layer are 

v = §3*5, (1) 
^ + (2-q)u = -d y Z, (2) 

^ - 2SO z , (3) 
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(4) 



where x, y represent the non-dimensional radial and az- 
imuthal coordinates and t is the non-dimensionalized time 
(Umurhan, 2008). In this construction, all quantities ap- 
pearing are functions of these two coordinates and time. 
The radial velocity is u. The total time derivative is 

— -> d t + (v - qx)d y + ud x . 

The background shear = —qxy, is linear in the radial coor- 
dinate. For a Keplerian disk, q = 3/2. Thus, the azimuthal 
velocity v appearing here is understood to be a deviation 
from this basic background state as viewed in the local ro- 
tating frame. The quantity H may be understood to be the 
enthalpy of the disk midplane and is given by 



E=ih 2 , 



(5) 



where h = h(x, y, t) corresponds to the height of the disk 
measured from the disk midplane. The height h is the po- 
sition at which the enthalpy goes to zero in the constant 
(specific) entropy environment. Although 2 is technically 
the midplane enthalpy, it is referred to hereafter as sim- 
ply the enthalpy with full understanding that it is really its 
value at the midplane. 

Equation (fTJ) is the radial momentum balance equa- 
tion but in this asymptotic limit it says that all pertur- 
bations are radially geostrophic. Equation ([2} is the az- 
imuthal momentum balance and, unlike the radial version, 
no geostrophic balance is implied. Equations (|5J) and flU 
represent, respectively, the motion of the surface h and an 
equation of state for the gas (see below) . In the derivation 
contained in Umurhan (2008), these two equations result 
from (i) exploiting the hydrostasy of perturbations at all 
times and (ii) an analysis of the energy equation, which 
makes use of the linear independence of powers of z. 

The disk vertical velocity, though not explicitly appear- 
ing, is odd symmetric with respect to the disk midplane 
because of the assumed varicose symmetry of the disk re- 
sponse. In the development of the DSW-theory its func- 
tional form is shown to be w = zfi z (x, y, t). In the reduction 
leading to these equations only d z w appears as well as the 
explicit evaluation of the velocity of the layer at the height 
z = h. Additional details are found in Umurhan (2008). 

The combination of these equations of motion leads to 
the conservation of potential vorticity Q given by 



DQ 
Dt 



-7^7- = 0, Q 



2(2 -q) + d 2 x E 



1 7 + 1 
7 2 7-1 



(6) 



The next section considers the fate of normal modes for 
potential vorticity (or ! PV for short) profiles in which there 
is either a single or a double defect in the steady PV profile 
Q = Q{x). 

3. Analytical study 

The analysis in this work considers the fate of disturbances 
where 7 = 3. This is chosen because the resulting equations, 
both steady and normal mode disturances, may be analyzed 
without intensive use of numerical methods. General more 
realistic values of 7 require numerical evaluation of both the 
steady and time-depenedent state and this will be reserved 



iui a future study. For 7 = 3, the PV expression takes on 
the quasi-linear form 



for 



2(2 -g) + 0^5 



(7) 



The following sections consider the two analytically 
tractable defect profiles and their normal mode responses. 

3.1. Single defect 

In this study, the term defect is used to indicate places 
where steps occur in the mean PV-profile. The following 
simple single defect, 



1, x<0, 
Qo x > 0, 



(8) 



is studied. The value of Qo is constrained to be greater than 
zero. Solutions are sought that are bounded as x — ¥ ±00. 
This means solving for the mean enthalpy S 



-1. 



(9) 



on either side of x = (and where q = 3/2 has been ex- 
plicitly used). The particular boundary conditions for the 
profiles at infinity are (i) S — > 1 as x — > — 00 and (ii) 
S — ► 1/Qo as 1 -> 00. At the interface x_= 0, it is also 
required that 3 and its radial gradient d x E are matched, 
the latter being identical to requiring the matching of the 
azimuthal velocities v from both sides. The solution for the 
mean enthalphy is therefore 



1 + O e x , x<0, 

1/Qo - e /VQ e-v^o* x > 0, 



where 



e 



l-Qo 
Qo 



(10) 



(11) 



This result may be read to mean that the mean height h 
decreases when Q Q > 1 and increases when < Qo < 1- It 
follows from Eq. Q that 



Ooe x 



x < 0, 



e oe -v^ - ] x > 0, 



(12) 



which means that the deviation velocity at x = is given by 
vo = Bo/2. As x — ¥ ±00, the mean velocity deviation goes 
to zero. The mean states developed here are qualitatively 
sketched in Fig. [1] 

Infinitesimal perturbations of this steady state are intro- 
duced by writing for each dependent quantity F — > F + F' , 
where F'(x,y,t) is small compared to the basic state F. 
This means that the potential vorticity equation is 



d t + (v - 3x/2)d y Q' + u'^- = 0, (13) 



where Q is given as 
Q> 

Because dQ/dx entails a delta- function at x = 0, spe- 
cial well-known treatments must be utilized to solve these 
equations (for example, see Drazin & Reid 1982). For this 
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type of problem, the procedure is detailed, for example, in 
Umurhan (2008), and may be summarized as follows: (i) 
solve Q' = on either side of the defect, (ii) match 5' com- 
ing in from both sides, and (iii) make sure that the normal 
velocities vf are matched across the boundary. The last of 
these may be straightforwardly enforced by inspecting the 
linearized version of Eq. @ , 



(d t + (v- qx)d y S j d x E' + (4 - 2q + d*E 



-23 



(14) 



and requiring that the two solutions coming in from either 
side of x = satisfy the requirement that the normal ve- 
locities u' are the same there. 

Normal mode solutions of the form H' = 
5(x) exp ia(y — ct) are assumed. Solutions with c > 
correspond to waves that propagate in the forward y 
direction and, correspondingly, in the negative y direction 
for c < 0. It follows that the solution for E/, away from 
the two boundaries, satisfying the condition that their 
amplitudes are equal at x = and, S — > as x — > ±00, is 



S = A 



e x , x < 0, 
x > 0, 



(15) 



where Aq is an arbitrary constant. The requirement that 
the normal perturbation velocities, u' , match at the defect 
location means ensuring that 



{vq-c- 3x/2)d x ~ + 2E 
1 + 81E 



(16) 



Satisfaction of this condition leads to an eigenvalue condi- 
tion on the wavespeed c which, after some algebra, reduces 
to the simple form 



v = -26 = -4u . 



(17) 



For vq given above, it follows that that c = — 3vq. The 
wave here is a Rossby edgewave phenomenon (e.g. Baines 
Mitsudera 1994) and the eigenvalue determined here de- 
scribes the speed with which the edgewave propagates. If 
the defect were placed at some position other than x = 0, 
then the observed speed of the edgewave would be Doppler- 
shiftcd according to the Keplerian velocity at the position 
of the defect. Thus, if the defect were instead at x — xo, 
the wavespeed would then be given by c = —3xq/2 — 3vq- 




a - Separation Scale 



Fig. 2. Plot of £, the ratio of the minimum or maximum of 
the mean enthalpy to its value as x — > ±00 (symmetric step 
profile only). Three values of Qo are shown: Qo = 0.25, 0.5 
(depression) and Qo = 3.25 (peak). 



half-width. Following the same procedures as outlined in 
the previous section and requiring that 3 — > 1 as x ±00, 
the mean enthalpy is found to be 



6 e^ 



1, 



- = < 



Bo cosh \/Qgx 
VQo sinh \/C7o A ' 



e e 



1. 



x < -A, 
-A < x < A, 
x > A. 



(19) 



The value of 0o in this case is 

1 — Qo ^/QotanhVQoA 



e 



Qo 1 + VQq tanh V'QoA 



(20) 



It means that for < Qo < 1, the enthalpy deviation is a 
peak or bump, while for Qo > 1 the enthalpy structure is 
a depression. The corresponding mean velocity deviations 
are 



6 e 



x+A 



x < -A, 



"=2 < 



-e ^#f, -A < x < A, 



(21) 



-e oe 



-x+A 



x > A. 



3.2. Double defect: The symmetric step 

This subsection is an analysis of a symmetrical defect pair. 
The previous section demonstrates how a single defect pro- 
duces a single Rossby edgewave along the defect. For two 
defects, therefore, two Rossby edgewaves are expected and 
it is under these circumstances that instability may emerge 
according to the principles of interacting waves (Hayashi & 
Young 1985; Baines & Mitsudera 1994; Heifetz et al. 1999). 

The mean PV-profile appropriate for a double symmet- 
ric defect is 

f 1, x<-A, 
Q = < Qo, -A < x < A, (18) 
[l, x>A, 

where A > 0. The edges of the defect are located at x = ±A 
and separated by a distance 2A, where A is the region's 



Defining w(±A) = v(x — ±A), it follows that at the defects 
v(±A) = =fOo/2. It is useful to define £, i.e., the ratio of the 
height of the peak mean enthalpy to its value as x — > ±00, 

max(S) = J_ _ _8q 1 

S(±oo) Qo vQosinhvQoA 

A qualitative sketch of the mean states are depicted in Fig. 

m 

Because there are two defects, individual edgewaves 
propagate along each of them, which interact with each 
other as well. A normal-mode analysis, as executed in the 
previous section, is performed. The details are presented 
in the Appendix. The resulting wavespeed c satisfies the 
relationship 

c 2 = £, (23) 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

width: A 



Fig. 3. Stability boundaries for the double defect problem 
in the A — Q plane. Unstable region depicted by shaded 
zone. Instability occurs for all values of enthalpy bumps 
(i.e., PV-depressions: < Qo < !)• No instability expected 
for enthalpy depressions. The boundary A = A cr (Qo) sep- 
arates stable and unstable regions. Critical values of Qo 
corresponding to instability goes like A -2 for A ^> 1. 



where 




+ 



1-Qo (cosh/3 + VQosinh/3) 2 

4 (/3 + V^o) sinh/3 [(Q„ - 2) cosh/? - VQ^sinh^] 



2Q cosh 2/3 



<3o)sinh2/3 



and j3 = y/QoA is defined for convenience. The results of 
this are plotted in Fig. [3] The first clear result is that 
instability (£ = 0) may only occur for enthalpy bumps, 
which is to say that instability occurs for PV-depressions, 
< Qo < 1- The boundary between stable and unsta- 
ble parameter regimes as a function of Qo is given by 
A cr (Qo). Thus, for PV-depressions, instability occurs for 
< A < A cr . For very small but non-zero values of Qo, 

— 1/2 

the results show that A cr ~ Q . Defining instead Q cr 
as the marginally stable value of the PV-depression then 
its asymptotic form has the functional dependence ~ A~ 2 
for A 3> 1 . Figure H] shows plots of the wavespeeds as a 
function of A for representative values of a PV depression 
and peak. 



4. Comments and discussion 

It has been shown that the RWI is recovered in the frame- 
work of DSW-theory. The RWI emerges when the idealized 
symmetric-step mean (radial) potential- vorticity profile is a 
depression (or equivalently, when the mean radial enthalpy 
profile is a peak). However, single stepped mean PV pro- 
files do not support the instability. The RWI is interpreted 
as being caused by the interaction of two Rossby-edgewaves 



Fig. 4. wavespeeds as a function of A for symmetric step 
profile: (a) Qo = 2.5, (b) Qo = 0.5. Also shown are the 
speeds c± (discussed in the text, Sect. 4) corresponding to 
the edgewaves that would propagate along each defect as 
if each defect were treated in isolation. Instability occurs 
in (b) only when these two wavespeeds are nearly equal, 
which is consistent with the Hayashi- Young criterion for 
wave instability. 



propagating azimuthally along each defect of the symmet- 
ric step PV-profile. These types of dynamical processes are 
a generic feature of many vortical dynamics in geophysical 
flows. Below is a collection of comments about these results 
including comparisons between the results of this work and 
the original RWI studies: 

1. A weakness of the analysis performed here is that an 
unrealistic value of the ratio of specific heats was used: a 
value of 7 = 3 places the fluid in a thermodynamic regime 
somewhere between a normal gas and an incompressible 
fluid (the latter being 7 — > 00, Salby, 1989). The main rea- 
son for this choice was to develop the analysis analytically 
as most other values of 7 require numerical methods to 
study both the steady states and the perturbation evolu- 
tion. Among other things, developing analytical solutions 
are important for benchmarking future analyses for more 
reasonable values of 7. But most of all, the analytical solu- 
tions offer the most transparency in interpreting the results. 
Preliminary calculations for realistic values of 7 show that 
the main mechanical interpretations and qualitative solu- 
tions, as well as conditions for onset, are not considerably 
altered. This is the subject of a study in preparation. 

2. The limiting form of the DSW-theory used to study 
the RWI considers the response of a single homentropic 
disk layer. In this picture, the mean vertical height of the 
disk is a function of the radial coordinate, i.e., h — h(x). 
Instability only occurs when the enthalpy of the disk layer 
(which is cx ft. 2 ) has a peaked structure. In LFLC-2000, the 
RWI is demonstrated in a vertically integrated model of a 
disk. Similar in quality to the results here, instability occurs 
when there is a peak in the mean radial pressure profile (or, 
relatedly, the surface density profile). In contrast, however, 
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instability occurs for constant specific entropy profiles. This 
might imply that these results conflict with one another. 

Because the disk model considered in LFLC-2000 is ver- 
tically integrated, all of the fluid quantities considered there 
are also understood to be vertically integrated, including 
the entropy. In the DSW-theory framework studied here, 
despite the specific entropy being constant throughout the 
layer, its vertical integration is not constant if the mean 
height of the disk is not uniform. That is to say if one de- 
fines 



-h 



Sdz = 2hC v \nK, 



as the vertically integrated specific entropy of the disk, then 
it is self-evident that S is not uniform if h varies with radius 
x. As such, the results of this study should be compared to 
the non-homentropic incarnations of the RWI discussed in 
the original studies, for example, as in Li et al. (2001). 

3. With respect to the last point, the version of the RWI 
uncovered in the analysis of the symmetric step is qualita- 
tively similar to the homentropic Gaussian bump profile 
('HGB') considered in LFLC-2000, where HGB describes 
the quality of the radial dependence of the mean pressure. 
However, as stated above, the symmetric step profile stud- 
ied in this work is effectively non-homentropic when viewed 
in terms of the vertically integrated disk model consid- 
ered in Li et al. (2000). The symmetric step profile might, 
therefore, most closely correspond to a non-homentropic 
Gaussian bump profile ('NGB') in the framework of the 
vertically integrated model of LFLC-2000. A NGB profile 
was examined in the nonlinear regime in Li et al. (2001), 
and, in a qualitative sense, its evolution was not found to 
differ considerably. 

One may, therefore, reasonably compare the conditions 
in which the RWI arises in the HGB setting to the instabil- 
ity arising here for the symmetric-step profile. In the analy- 
sis presented in LFLC-2000, it was shown that for the HGB 
profile a minimum amplitude in the pressure bump was re- 
quired for instability to occur. It was indicated in Fig. 9. of 
that study that the minimum variation required in the as- 
sociated surface mass density for instability to occur, on a 
radial scale similar to the thickness of the disk, must be be- 
tween approximately 10% to 20% (LFLC-2000). According 
to this study, for any given value of < Qo < 1, instabil- 
ity occurs for < A < A cr . This means that for a given 
peak value of the midplane enthalpy there is a maximum 
width size beyond which instability ceases. In this sense, the 
results here imply that there is no minimum required am- 
plitude of the midplane enthalpy peak (and, consequently, 
the corresponding associated surface density peak) for in- 
stability to occur. 

This does not necessarily mean that these results dis- 
agree with each other. Indeed in Fig. 3, it can be seen that 
for fixed values of A, there is a value of Qo smaller than 
some critical value for the system to be unstable. For exam- 
ple, for A = 0.7, instability occurs for Qo < 0.5. According 
to Fig. 2, this means that for these parameter values the 
minimum midplane enthalpy peak variation over the mean 
must be at least 25% for instability to occur. 

This still leaves open the question about the origin of 
the initial equilibrium. In the general scenario proposed in 
LFLC-2000, matter is accumulated in a piecemeal way, ei- 
ther by means of the coupling to large-scale magnetic field 



processes or due to small-scale turbulence (possibility MHD 
driven). These possibilities are not disputed here. Instead 
attention should be paid to the way in which matter is accu- 
mulated. Reference to Fig. 3 best illustrates this point. The 
way that matter accumulates depends upn the path within 
the Qo-A plane along which the accumulation process takes 
place. For example, if the process involves an increase in A 
(the region size) for fixed values of < Qo < 1, then insta- 
bility will occur before any sizeable enthalpy peak is reached 
because any such path (starting from A = 0) is unstable 
to the RWI. On the other hand, if the initial equilibrium 
Qo is constructed on a path of constant A, then the system 
may be brought into the unstable regime with a finite large 
amplitude of the enthalpy peak variation. Thus, the way in 
which the equilibrium is constructed is crucial to assessing 
how and on what scale the instability will operate. This 
point deserves further attention in the future. 

4. The single step profile in the mean PV considered 
in this analysis shows no instability. In contrast, the sin- 
gle jump profiles studied in LFLC-2000 do show the emer- 
gence of instability. The reason for this difference is not 
entirely clear. One possibility is that the single step profile 
considered here may not be a faithful analog of the single 
jump profile studied in LFLC-2000. Indeed, in the general- 
ized picture of edgewave dynamics of Baines & Mitsudcra 
(1994), one needs a minimum of two edgewaves for insta- 
bility to happen. The single-defect model of this study can 
only support a single edgewave and, therefore, no insta- 
bility is possible. A corresponding analog representation of 
the single jump profile considered in LFLC-2000 probably 
needs to support at least two edgewaves in order to see the 
instability expressed. This is, however, conjecture at this 
stage and more analysis is needed. 

5. Instability can be best rationalized in terms of the 
interaction of the two individual edgewaves propagating 
along each defect of the symmetric step profile. This pic- 
ture is consistent with the Hayashi- Young criterion for wave 
instability (Hayashi & Young 1987), which can be summa- 
rized as follows: instability of two waves separated by an 
evanescent layer occurs when the two waves (i) propagate 
opposite to each other, (ii) have almost the same Doppler- 
shiftcd frequency, and (iii) can interact with one another. 
Interaction in this case means that the edgewave at one de- 
fect contributes to the perturbation velocity profile at the 
other defect. 

All of these conditions for wave instability are satisfied 
in the problem encountered here. In particular, one can 
view the Rossby-edgewaves propagating along each defect 
as though they were disturbances of isolated defects. Thus, 
with respect to the wavespeed calculation worked out for 
single isolated steps in PV (Sect. 3.1), one can write down 
the wavespeeds associated with each defect as 



3 A 5 ^ 3 A 51 

c- = xA - -©o = -A - - — 
2 2 2 2 

for the wave at x = — A and 



3 A 5 ^ 3 A 51 

= — A + -6 = — A H 

2 2 2 2 Qo 



Qo 



/Qo' 



l 



/Qo' 



at x = A. The wavespeeds c± and the values of c appropri- 
ate for the double-defect problem are plotted in Fig. HI It 



3 The single wavespeed formula quoted above for c_ is written 
in terms of the formula given in Eq. (jlTf) augmented by the 
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can be seen in Fig.0t> that instability occurs when the pos- 
sibility exists for the two isolated waves to have nearly the 
same (and opposite) wavespeeds. In the rest frame x = 0, 
this occurs when each of the wavespeeds are nearly equal 
to zero. In Fig. 0k., where instability does not occur, the 
wavespeeds are not nearly equal enough (they are not close 
to zero), hence instability cannot occur. 

In this particular double-defect setup it follows that one 
can roughly predict the onset of instability to be when the 
isolated wavespeeds are nearly zero in the reference frame 
of an observer at x = 0. Using the formula for c± given 
above, c± ps happens very nearly when 

^ ^ 5 1 — Qo \/Qo 
~ 3 Q ' l + VQ^' 

Clearly this can never happen when Qo > 1 as A must 
be positive. Furthermore, when Qq <C 1 this rough crite- 
rion means that the critical value of the separation scale is 
approximately 

A ~0~ 1/2 

which is the same trend observed in Fig. [3J 

The importance of this perspective is that one may ap- 
proximately assess whether or not a more complicated look- 
ing PV-profile can become unstable to the RWI instabil- 
ity by appealing to this conceptually simple mechanistic 
view. For example, if a more sophisticated PV-profile were 
treated as if it were a series of PV-steps, then an approxi- 
mate prediction for the onset of instability would only re- 
quire calculating edgewave speeds along each defect taken 
as though each were in isolation. Onset conditions would 
be predicted for values of the parameters in which any two 
of these wavespeeds are nearly equal. 

6. The asymptotic limit that the DSW-theory exploits 
is one in which the azimuthal length scales are much longer 
than the radial ones. In this limit there is no information as 
to how wavespeeds depend upon the horizontal wavenum- 
ber a. In the theory developed here the wavespeed c is in- 
dependent of a. However, it is known that in the vertically 
integrated configurations investigated in the original RWI 
studies there is a clear dependence of c with a. 

On more concrete terms, in the numerical solution 
shown in LFLC-2000 there is obviously a fastest growing az- 
imuthal wavenumber of the instability, which leads to the fi- 
nal steady pattern state reached by the system. In contrast, 
in the approximate study done here for the symmetric- 
step profile, once stability is breached all wavenumbers be- 
come unstable and there is no high-wavenumber cut off. 
The growth rates, a = alm(c), obviously increase without 
bound in proportion to a - signaling the breakdown of the 
theory. 

To correct this issue, the DSW-theory would have to be 
appropriately extended by allowing for vertical structure in 
the horizontal velocities (i.e., inclusion of the dependence 
~ z 2 ) and some deviations from the linear dependenc of the 
vertical velocity with the z coordinate (i.e., an additional 
~ z 3 dependency). The resulting next order correction to 

Doppler-shift consistent with the the local Keplerian speed at 
x — —A. The isolated edge wavespeed c+ is correspondingly 
written taking into consideration the symmetry of the mean 
PV-profile as viewed at x = A, which is the opposite of that at 
x = -A. 



the wavespeed in the extended DSW-theory will show that 
the general dispersion relationship takes the form 

where £^ is a function of the parameters Qo and A. 
Preliminary calculations show that £^ is greater than zero, 
which means there is a high-wavenumber cutoff beyond 
which there would be no instability in the extended the- 
ory. The DSW-theory may be extended to higher orders in 
e by asymptotically developing corrections for c near val- 
ues of £ ps 0, which occur for values of A near A cr . Further 
investigation is needed in this regard. 

7. Three-dimensional numerical simulations of the RWI 
were announced during the preparation of this manuscript. 
Meheut et al. (2010) report on three-dimensional simula- 
tions of a constant entropy disk with midplane symmetric 
perturbations with the main goal of following the evolu- 
tion of a Gaussian bump profile in the initial pressure field. 
Notable features of the results show that there is a sig- 
nificant amount of three-dimensional structure in the hor- 
izontal velocities as well as deviations from a linear z de- 
pendence of the vertical velocities. The resultant structures 
show interesting and complex three-dimensional structure 
in the streamlines of the flow. In contrast, in the asymptotic 
three-dimensional theory used here the horizontal velocities 
have no vertical variation and the vertical velocity is linear 
in z. Thus, the asymptotic strategy used here is far from 
matching the overall structure contained in the results of 
these particular simulations. The DSW-theory (or an anal- 
ogous approach) may be extended to capture the quality 
of the dynamics contained in these numerical simulations 
much in line with the previous comments made in this 
regard. Moreover, it is important to understand how the 
three-dimensional structure emerges in these simulations, 
i.e., whether it arises from other secondary barotropic pro- 
cesses or if they are dynamics that are tied to the primary 
quasi-two dimensional dynamics driven by the primary in- 
stability. Rationalization of the three-dimensional dynamics 
in terms of the Rossby edgewave perspective utilized in this 
work could be a natural first step to this end. 
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Appendix A: Double-defect normal-mode 
calculation 

Following the procedure performed for the single defect, 
the equation Q' = is solved separately in the three zones 
subject to the condition that S — > as x — > ±oo. For 
x < — A and x > A, the equation to be solved is 

d 2 x 3 - 5 = 0, (A.l) 

while for —A < x < A, it is 

d 2 x E - Q E = 0. (A.2) 

The solution for S in each of the three zones is 
A_e x+A + 1, x<-A, 



§(A + +A_) co t^gs" 

2 v T ' cosh vQo A ' 

- A_) si1 ^^ , -A < x < A, 
A +e ~ x+A + l, x>A, 



(A.3) 



The solution as presented satisfies the continuity of S at 
the boundaries x = ±A. The remaining conditions are the 
continuity of the normal velocities at the location of the 
defects, which amounts to simultaneously satisfying 



(«(A)-3A/2-c)a x H + 2S 



1 + ^5 
and 

(w(-A) + 3A/2 - c)d x E + 2E 



dlE 



= 0, (A.4) 

A- 



= 0. (A.5) 

>-A- 



These two equations result in the matrix equation 
MA = 0, 

where M is a 2 x 2 matrix and A = A-) . Nontrivial 
satisfaction of this matrix equation requires enforcing 
det M = 0. This results in the wavespeed relationship found 
in Eq. (gl. 



